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SUMMARY 


The method of characteristics has "been linearized hy assuming that 
the flow field can he represented as a basic flow field determined hy 
nonllnearlzed methods and a linearized superposed flow field that 
accounts for small changes of boundary conditions. The method has been 
applied to two-dimensional rotational flow "where the basic flow is 
potential flow and to axially symmetric problems "where conical flows 
have been used as the basic flows. In both cases the method allows the 
deteiminatlon of the flow field to be slD 5 )lified and "the numerical work 
to be reduced to a few calculations. The calculation of axially sym- 
metric flow can be slnplifled if tabulated "vnlues of some coefficients 
of the conical flow are obtained. The method has also been applied to 
slender bodies without symmetry and to some three-dimensional wing 
problems "where two-dimensional flow can be used as the basic flow. Both 
problems "were unsolved before in the approximation of nonllnearlzed flow. 


INTRODUCTION 


The use of the method of characteristics for the solution of 
supersonic-flow problems requires numerical procedures "which are lengthy 
and Involved and "which must be repeated for each set of boundary condi- 
tions. The method has received general practical application only for 
two-dimensional or axially symmetrical problems in steady flow and dne- 
dimensionaJL or quasi- one-dimensional nonsteady flow, and only very few 
cases of general three-dimensional flow have yet been investigated. 

Many problems have been investigated at present by means of the 
linearized "theory in "which the disturbance-velocity con^jonents (defined 
as the difference bet"ween the local and "the free-stream con 5 >onents of 
the velocity) are considered and are small, so that terms of second 
order or hl^er can be neglected. In the present paper a sinq)llfication 
is introduced in the equations of motion based on "the assuu 5 )tion that 
one of the veloci"ty con^jonents or the variation of the velocity com- 
ponents as a function of a given parameter can be considered small, so 
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that terms of second or higher order in the quantities considered small 
can he neglected. When one of the velocity components is assumed to he 
small, the other two velocity con^jonents can he e 35 pressed in two parts, 
one of which is large and is a function only of two coordinate positions, 
and the other of -vdiich is small, of the same order as the third velocity 
con5)onent, and is a function of all three coordinates. 

If the variations of velocity conponents as functions of a given 
parameter are considered small, all three velocity conponents can he 
estpressed in two parts. One, large, is independent of the para m eter 
considered, and the other, small, is a function of the parameter con- 
sidered. When the velocity components are substituted into the differ- 
ential equations, the equations can he divided into two parts, and the 
differential equations containing the velocity components considered small 
hecome linear; therefore, supeipositlon of solutions is possible. With 
this assumption the flow field can he represented for any condition by 
the superposition on a nonlinear basic flow field of a linearized flow 
perturbation. The flow field, 'vdiich represents the variation of the 
basic flow due to the changes of the geomietrical parameter considered, 
changes linearly with the parameter. Because of the simplification, 
the superposed flow field is defined by differential equations of hyper- 
bolic type which have characteristic surfaces equal to the characteristic 
surfaces of the basic flow field and Icnown coefficient; therefore, the 
supeiposed flow field can be obtained directly without the iteration 
process along the characteristic net of the basic flow. 

A particular application of the linearized characteristics method 
has been discussed in references 1 , 2 , and 3 in which bodies of revo- 
lution at small angles of attack have been considered. In the present 
paper the basic concept of the linearization is discussed and examples 
of application to two-dimensional rotational flow, to conical flow, to 
axially symmetric flow, and to some general three-dimensional problems, 
are discussed. From these examples, other applications of the same method 
to supersonic steady- or nonsteady-flow problemis can be visualized. For 
example, the miethod can also be applied to the determination of the flow 
field in supersonic compressors or turbines having supersonic relative 
velocity inside the passage. In this case, the two-dimensional flow of 
the cascade of the blades at each radial station or the axially sym- 
mietric flow can be assumed as the basic flow. In the first case the 
radial component of the velocity must be assumed to be small, while in 
the second case, the tangential component of the velocity must be 
assumed to be small. 


SYMBOLS 


‘a,b,c,d,e 


coefficients of the variables in the characteristic equa- 
tions (defined case by case) 
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a 


speed of sound 


An 

B and C 


M 


coefficient 

coefficients defined by eq^uations (20) and (25) 
pressure coefficient 
Mach number 


R 


gas constant or reidius of hodograph diagram 


S 

u,v,w 




entropy 

velocity components along the x-, y-, and z-axes in Cartesian 
coordinates or along x- sind y-axes and peip)endicular to 
the meridian xy-pleine in cylindrical coordinates 

velocity con^onents in polar coordinates 


V 


tangential and normal velocity components in front of the 
shock 

Intensity of the velocity vector 




Mach angle 


7 


ratio of specific heats 




inclination of the characteristic lines in the plane 
z = Constant or 0 = Constant 


0 


coordinate of the meridian xy-plane in cylindrical coordi- 
nates or of the rt-plane in polar coordinates 


^,9 


polar coordinates 




inclination of the velocity vector with respect to the 
x-axis 


components of the rotation along the x-, y-, and z-axes 

Subscripts : 


o properties of the basic flow 

l...n properties of the supeip)osed linearized flow fields 
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THE EQUATIONS OF THE LINEAEIZED CHAEACTERISTIC SYSTEM 


Consider, for exan^jle, a flow field defined ty a velocity vector 
V(u,v,w) the con5>onents of which can he expressed in the form 



( 1 ) 


w 





The relation between entropy and rotation states 

— — 2 

curl V X V = — grad S 
7R 


( 2 ) 


Assume that Uq eind v^ are functions only of x and y, then hy 
neglecting terms of the order of a^^, eqization ( 2 ) becomes 


Sx 7R °\Sx by } 




^ = u 

'by 7 R 




^ ( 3 ) 
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Tlierefore, 


the entropy field can be expressed as 


1 




> W 


Sz 



bSn 

bz 


Assume that each coefficient a^ . . . is constant in the 
entire region of the flow field where it is not zero and is small, so 
that terms of the order of a^^ or higher can be neglected. In this 

flow a basic flow field exists, represented by the velocity 
vector Vq(uq,Vq^ end by the entropy distribution Sq, on which a 

linearized flow is superposed, represented by a summation of N three- 
dimensional flow fields, each of ■vdiich is proportional to the corre- 
sponding coefficient a^. The basic flow field is a two-dimensional 

flow if Cartesian coordinates are used or an axially symmetric flow if 
cylindrical coordinates are used and can be determined by known methods. 
In a similar way, a general three-dimensional flow field can be assumed 
for basic flow if the flow field can be obtained by sln5)le analysis. 

The equation of motion obtained from continuity, momentum, and 
energy equations can be expressed in Cartesian coordinates in the form 


^4 _ u£) + ^ I _ 

5x \ a2y Sy\ a^j ^z\ a^J a^y^y ^xj 

uv/bu. ^ 5w \ _ wfbv _ Q 

a^ySz bx) a2\^z by) 


( 5 ) 
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while, in cylindrical coordinates, the equation hecomes 


^(1 - + ^(1 - + ^2-/1 uv /Su ^ ^v\ 




By\ a^y y a^/ a^\Sy Sx 


wu | 5u . 5w \ _ wy|_^T_ . + X. = 0 

a^ly Bcp Sxl a2W Sep Syi y 


(6) 


By use of equations (l), expression ( 5 ) hecomes 


^(1 _ liofV _ 

V a 2] Sy i-" .2 


Ur.V, 


Qlq/XXo 


ao^/ ao 2 \Sy Sx 


Ur^ 5V, 


z 


«n 


^ + ^/l _ 


hv,. 


Sz Sx 


+ _ 10, 




lio^o/za + ^ 

ao 2 Sy J 


aup /2upUn + Up2Ap\ _ 

I a„2 i ■ Sy \ 


’-o 


/Suo 

'^o^cAi + ^o^n + ^ 0 % 

\Sy ^x J 

a 2 

ao 


( 7 ) 


where 


An = + ^o-^n) 

®-o 


The basic flow is a two-dimensional flow, and the first part of 
equation (7) must be zero; therefore, the second part of the equation 
must also be equal to zero and, for each value of n. 
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^ 

Sz Sx \ a.Q^j hy \ slo^J &q^ \Sx hy J 
Buo/2uoUn + Uo%^n\ &Vo|^oVn + Vo^AjA 

;%y\ ao^ r 

+ ^^oVo^cAn + ap^n + ^0^11 
V^y a^2 


Equation (8) is a differential equation of hyperbolic type because 
Uq^ + Vq 2 > and the characteristic surfaces are cylindrical surfaces 

perpendicular to the plane z = Constant with generatrices coincident 
with the characteristic lines of the t-wo-dlmenslonal flow. This can be 

seen from the fact that the coefficients of the derivatives 

ox dy 

8v„ 8u_ , - . 

and — ^ of equation (8) are expressed as functions of the prop- 

ox 8y 

erties of the basic flow and are the same as the coefficients of the 


derivatives 


and which define the basic flow. 


8x " ^y , ^ 

8v_ 

and the coefficient of g-= is one (see, for exanple, reference 

page 282). Therefo3?e, the disturbance flow field 7 jj(u^,Vjj^,Wq) can be 

obtained by the method of characteristics by moving along the charac- 
teristics surfaces which are cylindrical. The characteristic net, which 
for the general case must be drawn in spatial coordinates, can be drawn 
in this case only once for any value of a and n and is equal to the 
net of the two-dimensional flow. Equation (8) can be transformed for 


practical use 
8vo 


In the equation, the terms in UoVq, ^Z2 

ox dy dx 


and 

dy 

field. 


are known terms and are given from the two-dimensional flow 


Assume, in the plane xy, polar coD5)onentB for V and cp for the 
velocity, defined by 


Uo = Vq cob cpo 
■^o = "^o 'Po 


(9) 
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ajid, ‘because M^er-order terms are neglected^ 


■wiiere 


Uq = Vji cos cPq - cPjjVq sin cp^ 
sin <Po + ^Pa^o % 


( 10 ) 


cp = cPo + a^cpj, 


J 


( 11 ) 


If n Is the normal to the projection of the streamline in the 

plane z = Constant, the derivative ^ is e:^)resaed hy 

on 


8S _ 8S V j 8S u 

8 n 8 x V 8 y V 


( 12 ) 


By use of equations ( 4 )'€ind (l). 


equation ( 12 ) becomes 


8 n 




By considering only the lowest-order terms, equation (13) becomes 


8S ^^o \ — 









( 14 ) 
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But, 


3 x 


L/'v V - T V ^ /'ll Y - u V = -cp V ^ Q 

o Wnj + Wo UoVnj Wo Vq 6 y Vq/ 


■because the ri^t term contains as a factor the variation of entropy 
along the streamline of the "basic flow; therefore, 


5Zn ^ = [Yr+ ^ 

ax “ ay rW ano \Y^ ^jr^o 


(15) 


where ng is the normal to the stream! in e of the "basic flow in the 
plane xy. 


After several transformations the following eq[uations can he 
obtained (see reference 2 ): In the plane z = Constant along the char- 

acteristic line defined "by 


dx 


Xg = tan(P^ + cp^) 


(16) 


the following equation is valid: 

d•^^ 1 sin pQ tan Pq ^ 1 dV^ ^ ^ sin^Pp ^ dS^ ^ 

dz Tp cos(Pp + cpp) Vp dx ° dx cos(cpp + Pp) 7R dnp 

'PnBl + ^ Cl = 0 (17) 

'' o 

and along the characteristic line of the second family defined hy 


g = X2 = 



(18) 
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the following equation is valid: 


dv^ 1 sin Pq tan Pq 1 , "^-Vn 1 , ton p 1 

dz Vq cos(cpo - ^Q) Vq dx Vq ° dx cos(cpo - 7K 

% 

+^<=2 = 0 ( 19 ) 

* o 

Tdiere and Cg are coefficient functions of x, y, 

cpQ, and along the diaracteristic lines of the first 

and second families at each point and are independent of q)^, 

and Wjj^, and, therefore, can he determ in ed once for the basic flow and 
used for any kind of disturbance flow in the limits of the approximation 
accepted. The coefficients B^^, B 2 > C-j^, and Cg are given by the 

following expressions: 


Bl 


1 

sinVo <iSo 1 


cos(cPo - Po) 

cos pQ cos(cpQ + Po) 

7R dUo Vq 

(dx j 

'2 sin Po 


(20a) 


Bg = 


cos Pq cos(cPo - Pq) 


sin^o dS^ 1 /dVo\ cos(q)p + p^) 

7R dn^ ^1 Po 


(20b) 


_ iK 


Cl = /tan% + 


7-1 


^o\^ >± \ 2 sin^P^ cos^p. 


2 ^0 sin3p^ 3^1 

^dVo) 

cos(cPo - Po) 


^1 1 ^ ~ ^ \ 

7R too cos(cp^ + P^) Vo' 

^dx 

2 cos(cPo + Po)cos2p^ 

V 2 sin^p^y 


(20c) 
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C 




+ 


7 - 

2 sin% 



+ 



ala\ 


cos(cPo + Po) 1 

/i I ^ \ 

7R dn^ 

cos(cPo - Po) 

VoW /i 

cos(cpo - Po) cob \ 

( 2 Sin2p^) 


( 20 d) 


viiere 



is the derivative along the characteristic of the first 


family and 



along the characteristic of the second family. 


In order to determine the value of at each point of the char- 

acteristic net^ the following relations can he used in the plane 
z = Constant; 



bv 

Sy 

bz 

du _ 

bv 

Bz 

dx 

bv _ 

Su 

Bx 



( 21 ) 


and, if Sq is a streamline projection in the plane z = Constant, 

Sw_ _^u,^^v^_u v 1 + a£_ 

3sq Sx V Sy V V V hz bz Sz 7 RV 


( 22 ) 
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In the approximation accepted. 


^ ^ 
Bz ^ ^ ^z 


Sz 



(23) 


therefore. 


_L ^o^ 

^o Vq Sz Sz yRv^2 


m 


•vdiere is the streamline of the basic flow. 


Equation (24) permits the determ in ation of the value of w at each 
point of the plane z = Constant as a function of the local variation 
of Vq in the z-dlrection. Therefore,' hy means of a step-hy-step pro- 
cedure, all the flow field can he determined hy working in planes 
z = Constant ’where only one characteristic net must he used, hy means 

of eq.uations (1?) and ( 19 ) and hy determining the value of at each 

dso 

point of the net hy means of equation (24). The calculations are started 
from the flovr field defined, or along a surface -which is not a charac- 
teristic surface, or along a characteristic surface and a stream surface. 
The flow field at the starting surface must he determined from the 
boundary conditions. If the starting surface is a shock -wave, the flow 
at the shock surface must he oh-tained from the physical properties of 
the shock -wave related to the boundary conditions considered. Relations 
het-ween boundary conditions and shock -waves are presented in detail 
subsequently for the problems considered. 


Similar' equations can he obtained hy using cylindrical coordinates 
in place of Cartesian coordinates. 


Equations (l?) and (19) remain the same. Only the first term 


Sz 


5w 

T 


and the coefficients and B 2 become 


changes. 


becomes 
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1 sin Po sln(cpo - p^) 


cos Po cos(cpo + Po) y 


®2* = 


sin po sln(cpo + p^) 


cos Pq cos(qp^ - p^) y 




Bg 


(25) 


Equation (25) for cylindrical coordinates ‘becomes (see reference 2, 
equation (9c)) 




I 

0=Constant 


y de 7B7^ 


(26) 


The use of the linearized characteristic system is sin^jler in many 
cases than the con^plete characteristic system and reduces in some prac- 
tical applications the extent of the numerical work required, especially 
if solutions of many similar pro’blems are req,uired. The same concept 
can also he easily applied in the field of nonsteady or relative motion 
of flows. 

The flow field around slender bodies without axial symmetry can he 
ohta±ued hy mjans of linearized methods, and some practical three- 
dimensional problems not a na lyzed before in the approximation of non- 
linearlzed flow can be analyzed by this method. 

In the next sections some typical possible applications are 
presented. 


SOME TYPICAL APPLICATIONS OF THE LIHEAKIZED CHAEACTEEISTICS METHOD 
Two-Dimensional Rotational Flow Fields 


Two-d i m e nsional potential flow permits hodograph solutions, and, 
therefore, any kind of two-dimensional supersonic potential-flow solu- 
tion can "be obtained in the hodograph plane; the numerical solution in 
the physical plane then requires only the construction of a character- 
istic net in order to find the position in the physical plane of each 
point of the hodograph plane. The solutions of problems in which 
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■boundary conditions are given only along a streamline are very sin^jle 
in the approximation of potential flow "because in this case the velocity 
is constant along a family of characteristic lines (single-wave flow); 
however, similar calculations for rotational flow are much more involved, 
because a step-by-step procedure is required for the solution on the 
holograph plane as irell as for the construction of the characteristic 
net. By means of the linearized characteristics method a rotational 
flow field can be considered as a modification of a potential flow 
field, and the linearized superposed flow is the flow which takes into 
account the effect of the presence of shock waves and the effect of 
rotation in the flow. 

Consider, for example, a two-dimensional profile which produces a 
shock wave at the leading edge (fig. l). If the profile is curved, the 
shock is curved, and the flow behind the shock is rotational. Assume 
that the flow field behind the shock can be expressed as 


or 


u = Uo + ui 


^ = ^0 + ^1 1 


V = Vo + Yi 


= «Po + ^Pl 




( 27 ) 


and 


S = 


■\diere Uq and Vq are the velocity components defined by a potential 

dV 

flow field, which in this case is a single-wave flow, -3—2. along the 

characteristic lines of the first family is zero, and u^, v^, Euad 

represent the flow field that takes into account the reflections 
occurring at the shock and the effects of the entropy gradient. The 
flow field represented "by and can be immediately determined, 

and Vq and cp^ are constant along the characteristic lines of the 

first family, which are strai^t lines. The characteristic net of the 
flow VqCPq can be drawn in a short time. Then, along the characteristic 

of the first family. 


dVn i - 
''o 


sin^S 

tan 3 o dcpi + dSi + dx^ + ^ dx^ = 0 ( 28 ) 
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■where- 


and. 


Bi = 


COB ^Q COS 


(<Po + Po) 


cos(cPo - Pn) 
Bln Po 


(29) 


= J:. COs((Po - Pq) _i / 

\dx Vo cos(cpo + Po) co^2p^^ 


7 - 1 ' 
2 s±n%^ 


(30) 


and., along the characteristic of the second, family, 

, sln^ Vt dV / -V T \ 

i dV + tan p^ dxpi + ^1 + ^ V“^(tan P^ + 2 ] = ° 

° ''o Vo 2 sin% cob \) 


(31) 


Equations ( 29 ) and (3I) can he simplified hy introducing -the gra- 
dients (Vq)^ and ((pQ + Pp^^ along the x-axis defined as 


(«. ■ (£■) 


{% * »o), “ 


y=0 

^(q>o + Pq) 

Sx 


(32) 


Jy=o 


It can he seen from figure 1 that, at any point A of ordinate y. on 

A. 

the characteristic line a, crossing the axis at A^, the -\rari- 

/dV„\ 

atlon along the characteristic line of the second family at A 

is 


(&), 


“ (Vox)a 


1 - 


cot| 

(l>0 + So), 

•“•0 

1 

cot| 

('Po " Po)^ 

■^0 

Va 

(fo + So)^_ 

Ao 


sin^(cp^ + Pq) 


( 33 ) 
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because and cPq are tbe same at A and at Aq. Therefore, 


1 ^ a 

VoNdx 1 


(31^) 


•where a and b are constants along each characteristic line of the 
first family. Tbftn equation (28) along the characteristic line of the 
first family becomes 


dVn 


— i - tan 8_ dflK + 
Vo ° ^ 


7R 


d% - cp]_ 


cos((Po - Po) a dy ^ 

sin(cp^ + cos Pq sin Pq 1 - by 


Zl/l + y - 1 \ a dy cos((Po - P^) 1 

^o\ 2 sln^oj^ - sln(q)Q + p^) cos^p^ 

Therefore, if 


sin2(cpo + Pq) ’'^o 

and 

t Mx 

sln2((pQ + 3o) 






(35) 


(36) 


the equation along -the first characteristic line becomes 


^ - tan Pq dcp]^ + dS^ + 9 ^ I d log(l - by) 


1 + -Z.. - ^ -\tan Pq ^ - d log(l - by) = 0 


2 Bln^Po 


(37) 
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and, along the second characteristic line. 


d.Y-1 

— + tan p dcp. + 
V o ^1 

’'o 


sin^Q 

7R 


dS^ 


+ 


Vo 




+ —Ljl1 \ 

2 cob^^I 


= 0 


( 38 ) 


All the coefficients of equations (37) and. ( 38 ) are constant along 
characteristic lines of the first family and can he calculated at few 
points on the x-axis. 


The coefficients of equations (37) and ( 38 ) are independent of cp^ 
and ^therefore, with one calculation from a point A and C of the 

net of figure 1 the. values of and cp^ at a point B can he 

obtained directly without the necessity of an iteration process, 

including also tejmis of the order of (Ax)^ in each step. Indeed, if 
all the quantities -which are -variable along -fche characteristic lines are 
e3q)ressed in -the form 


and 


then 


f 8a\ ^ ^ ^ + o(ax^) 

\OX/g \OX/j^ 


(39) 


®B “ ®A - 



(40) 


Consider now equations (37) and ( 38 ). At point C equation (37) has the 
form 
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■while at point B it has •the form 



then, from equation ( 40 ), 
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and also 



+ 



2 


Ax = 



+ 



/ C 


Ax + 



Therefore, equation 
the order of fxp - 


( 37 ) 


along the first characteristic line^ terms of 
being neglected, can be written in the form 



and eq.uation (38) along the second characteristic line, terms of the 
order of (xg - being neglected, can be written in the form 







r - 1 

2 sin^3_ cos^f 


N - \) 


= 0 


( 43 ) 
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Because of the posslhlllty of considering directly in the calculations 

terms of the order of (Ax)^ also, large steps can he used in the 
cheiracterlstic net and the effect of entropy gradients can he easily 
evaluated, after the basic characteristic net and a few streamlines of 
the basic flow have been determined emd the coefficients of equa- 
tions ( 39 ) (4o) calculated at a few points on the axis. 

All the coefficients are constant along each characteristic line 
of the first familyj therefore, in going from B to D, only the terms 
containing S]_ and y must be changed and the calculations are simpli- 
fied to some extent with respect to the rotational- flow characteristic 
calculations . 


Axially Symmetric Flow Fields 

For axially symmetric flow, the equation of linearized character- 
istics becomes: 


"I dV-i dcpi sin dS-i , Vt 


along the characteristic line = tan^^Q + (Pq) 


T dT-t dJh sin%_ dSn i Y-, 


(^ 5 ) 


along the line X 2 = tan^cp^ - where ^2 

defined by equations ( 25 ) and (20). 


The Introduction of equations (il4) and (45) simplifies noticeably 
the numerical, calculations without affecting sensibly the precision of 
the results. The practical, use can be as follows: A basic body shape 

is determined first by means of characteristic calculations, and the 
characteristic net is then obtained. The basic calculations must be 
extended in a region in front of the shock wave determined with the 
usual procedure (for example, reference 3)^ in a region Inside the 
body as shown in figure 2 in order to determine all the flow field 
necessary. Thus, a reduced number of points of the characteristic net 
are chosen at which the supeiposed flow field for each boundary condi- 
tion different from the basic shape will be determined. The number of 
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points required depends on the magnitude of the supeip>osed flow field; 
however, the number is usually small, because for each step between two 
points A and B the disturbance velocity can be expressed in the form 



= (Vl) 

^ \dx ; 


Ax + 



(46) 


also included because, in the dlffer- 

the coefficients of the differential 
equations are independent of the solution and are knovn at both points. 
If the entropy term is neglected, by applying finite-difference 

methods the velocity con^jonents at a given point g5 of the character- 
istic net can be expressed from the values at two points f5 g4 in 
the foim: 


■vdiere the term 




is 


^dx^ 


ential equations (44) an<^ (45), 



(Vl) 


1 - n 

g4 1 + n 



m -f L 
1 + n 



m - L 
1 + n 


(47) 




1 - r 
f 5 1 + r 


+ 



P + <1 + 
1 + r 



<I - P 
1 + r 


m 


('Pi) 


g5 


1 

a - £ 

1 + n 1 + r 


(V, ) — 

\ l/gl;. 1 + n 


(Y,) 

' ^'f 5 1 + r 



m + L 
g 4 1 + n 



P + a 
f 5 1 + r 


(49) 


■vdiere 




and 


(^ 1 ). 


are functions of the bound- 


g4 ' -^'f5 

ary conditions considered, idille the coefficients L, m, n, p, q, 
and r are functions only of the basic flow field and, therefore, must 
be determined only once for any boundary condition considered. These 
coefficients are: 
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L = 


m = 


(Vq tan ^o) + (Vq tan 

g5 




- V) 




“■( 


■ '^g5 - ='g4) 


- {<= 1 )^ 


P = 


(7q tan + (7„ tan 


' {^g5 - ^f5 




r = 


(^g5 “ "'f?) 


i’=s)^ * (<=2)f5 


1 

5 : 


> ( 50 ) 


J 


d5-i , ^ » 

If the entropy terms in are considered., tvo more terms in 

and (S^^) must "be considered In equations (47) and (48) and two more 
terms in (S^) (®l) 1 ^ equation (49) . However, 'because these 

O ‘ 

terms are small, they can usually he neglected in practical calculations. 


Because all. the coefficients are determined only once,, the deter- 
mination of or (tPi) for each "boundary condition is slii^jle. 

g5 cr^ 


'g5 


Points at the boundary can he investigated hy means of equations (47) ^ 
(48), and (49). Points on the shock can he analyzed in a similar 
manner. A practical calculation c6iu he performed in the following way: 
The "basic flow field and the number of points in -vdiich the superposed 
flow field will he considered having been determined, the new boundary 
conditions (shape of the body) are placed in the characteristic net 
(fig. 3)* If A is the point where the "basic body departs from conical 
shape and OBC is the new boundaiy condition, the flow field between the 
surface of the conical body OB eind the conical shock OE are known from 
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cone calculations; therefore, the values of V and cp at each point Ba, 
la, 2a, 3a-, nnd Ea of the "basic characteristic net are determined, and 
"by difference the values of V]_ and can "be obtained. From al the 

flow at the point L of the boundary can be determined from equation (48) 
■where cp^ is known at L, and the coefficients at L can be determined by 

linear interpolation between bl and cl. From the values of and cp]_ 

at al and LI, the corresponding values at bl are interpolated. Then all 
the ■values for the line b can be ob^tained. For the detentii nation of "V"]_ 

and cpj^ at a point F on the shock wave, the equations of the shock wave 

and equation (4-7) are "valid. At the point F the "values of cPq and Y^ 

are known, and from the equations of the shock "wave the "value of — as 

8 <p 

function of cp can be determined. (The "value of (p fixes the de"viation 
across the shocdc "wave.) Therefore, 

(V)p = (Vo)^ ^ (Ti)^ = (v,)^ + (|)(?,)^ (51) 


Then, equation (47), applied be"trween the points 4b and F, gives 


“ (I), (^i)p 


' ^ ^ - (^1) 


m + L 


m - L 
F 1 + n 


(52) 


and the "value of can be determined. 

The "work requi3red in the cadculation of the flow field for the 
"basic body und the determination of the coefficients L, m, n, p, q, 
and r can be reduced to a minimum if conical bodies are assum^ as 
basic bodies for "the calculations, bec:ause in this case "the basic flow 
is available in tabulated "values (reference 5) "the coefficients L, 
m, n, p, q, and r are functions only of the polar coordinate ijr. 

For conical flow, the calculations can be performed in the 
following "way: From conical-flow calculations, the values of v^ and 

Vq as functions of i[r are known, "where Vj, is the radial con^ionent 
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and ttie normal con^jonent of the flow field in polar coordinates 

refeirred to the limiting velocity (fig. ^4-). The following e3cpressions 
for Vq, and cpQ can he determined from the cone calculations; 


= V ^ + V ^ 

'' ''n ^ ''r 

1=1^^ sin^p, 

» o 

q>o = <p + 7 


> (53) 


tan 7 = ^ 


Therefore, 6 + cp and op - • are known as functions of ilr. From 

^ O 'O o 

a point A on the conical body, the characteristic line of the first 
family AC and of the second family AD can he drawn. The line is defined 
hy the expressions 


__ _ 


J-p - 7^ tan(p - <p)^ + tan(p - cp)^ 
^ - ^A 2 


In order to construct a characteristic net that requires only a small 
eimount of calculation, the points of the net are chosen along straight 
lines from 0 so that the conical property of the flow can he utilized. 
The net can he constructed hy fixing the steps along the body. 

When point B is chosen along AO, the point E along AC is determined 
hy drawing BE parallel to AF. From E and F (along OE and AD) the lines 
GE parallel to LF and FG parallel to EH can he determined, and the 
point G can he obtained.. From G and H, point M can he obtained, and hy 
proceeding in a similar way, a].l the characteristic net can he deter- 
mined. The coefficients of equations (50) ajre the same' for each point N, 
F, and E along the same radius , The tabulated values are given only in 
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the region between the body and the shock wave; however, the calcu- 
lations can be extended by means of conical calculations. For exan^jle, 
from, the following equations (see, for example, reference k, p. 2 ^ 3 , 
and the following pages) 



Vj. + Vj^ cot tJc 


1 - 



(r - 1) (1 



and 


( 5 ^) 


“ ('"“Vo 


(''"'VcHAt ° 'S* - (® - V)„. ^ 


^ (55) 






(■sdiere R is the radius of curvature of the streamline in the hodo- 
graph plane); therefore, the characteristic net can be extended to the 
outside flow. For conical flow the coefficients B]_’, B2S C^, and C2 
can be determined as functions of ijr. For conical flow, coef- 
ficient B]_* of equation (25) becomes 


1 

Bin Pp Bin(cPo - p^,) 

cos cos(cp^ + p^) 

y 


2^ cos((p^ - P^) 

Vo dijr \dx/j^^ sin 
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and 1360001168 


C 


1 


Vo dt 


2q 

cos Po 


/sin'=^3o + 


7-1 
2 sin^ 



+ 


~ Pq) 1 A , -7 - 1 \ 

Vo dijr cos(cPo + 3o) coB^Pol^ 2 sin^p^j 


tut 


and 



sin t sin^cpQ + p^ - i]r) 
y cos(cPo + Po) 



sin t sin^t - (p^ + Po) 
y cos(cPo - Po) 


( 56 ) 


( 57 ) 


Tliarafore, 


a 

dijr 



1 

cos Pq 


Fo 


sin Pq sin((Po - Pp) 
sin t sln(cp^ + Pq “ ^) 


+ 



I 1 Sln(t - (Pq + Pq) 
Tjr sin Pq sin(cpQ + Pq - i{f) 


( 58 a) 
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and 



Bo V, 





1 

cos Po 



sin pp sln((Po + Pq) 
sin If sin(ilf - cp^ + p^) 


1 sin(q?o + Pp - t) 
\<3-'*lf /ijr sin Pq sin^ilf - cp^ + p^^ 


(58b) 




cos 




sin^Po 


+ 


7 - 1 > 

2 sin^Poy 


1 A ^ 7 - 1 - ^o Pq) 

cos^Pq^ 2 sin2pjsin(cp^ + P^ - iJr) 



1 

1 ^ 

sin(cPo + ^o - t) 

cos^Po 

^ 2 Sin^p^y 

sln(t - (p^ + p^) 
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and 


£!o = Zn Vr + ^ 

at “ Vo 2v 2 

^ cvn 

“ (7 + l)(l - Vq2) 


If equations (58) are used, the coefficients of equations (50) can he 
ejcpressed as follows: 


= t(^g5 - V) ^g5) 

^ ~ I^(^g5 “ ^g^)('^g^ ‘^g5) 

^ ^ ("^g5 ■ ^*^ 5 ) (^f5 ^g5) 

"" " ^g?) 




> (59) 


Because the coefficients a, h, c, and d are functions only 
of ■'}'■ and of the free-stream Mach number M, they need to be calculated 
only once for different "values of and gl"ven in "tabulated form, and, 

therefore, the calculation of any flow field for "vdilch the basic flow- 
can be considered a conical flow can be reduced to the solution of a 
few linear equations "wi"th known coefficients. 

In order to try the method, "the flow around an ogive, as shown in 
figure 5, has been determined at M = 3*Ol6, by the method of character- 
istics and by the method of linearized characteristics. In figure 5 the 
usual characteristic net is shown, "while in figure 6 the linearized 
characteristic net and "the basic body are presented. As a first basic 
body the cone chosen is the cone "tangent to the apex having a cone angle 
of 12.5°. The cone chosen is not the most convenient because the values 
of cp^ at the end of the ogive are large, and a better approximation 
would be ob"tained if a cone of smaller cone angle would be considered as 
the basic body. A cone some"what different from the ogive considered has 
been chosen in order to have some information on the approximation of 
the method for sensible "variations of the shape of the body from the 
basic body. 
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Downstream of the characteristic line AB, the velocity com- 
ponents V]_ and cp^ would become large because at the surface of the 

body the con5>onent cp is quite different from the cnmp nnpn t. cp^ of 

the basic flow; therefore^ the flow determination in the region down- 
stream of the line AB has been considered as a new problem, defined by 
the flow along AB and from the stream3jlne that represents the body 
shape. In this region a new basic flow has been considered. Again, a 
conical flow field has been assumed as the basic flow. From the values 
of cp and V at B and A, and ftrom the order of magnitude of cp down- 
stream along the body, a conic:al flow field that would give a small 
disturbance con^ponent in this region has been selected. The cone chosen 
for the second part is a 5° cone at M = 3.077. The cone is entirely 
contained within the body considered in the region used in 1iie c:alcu- 
latlons as shown in figure 6, and the most convenient region of the 
conical flow field is used for the calculations. 

In order to pass from one basic body shape to the other, the com- 
ponents and 9^ along the characteristic line AB for the second 

basic body must be determined. This operation can be performed by 
detennlning the value of Vq’ + and 9^* + cp^^’ at the points a, 

b, and n of figure 6 for the first basic body, by Intei^olatlng the 
values at the points a*, b*, and n* between characteristics 1 and 2 
along the characteristic of the other family and then determining the 
new value of and 9^" for the second basic body along AB from 

the eapressions 


= V - Vq" 

= 9 - q>o" 

and 

V = Vo* + Vi* 

9 = 9o» + 9^» 

where Vq” and 9 q" are the values of V and 9 for the second basic 
body at the position considered and V^’ and 9 qS "tbe values for the 
first basic body. 
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The calculations by linearized characteristics required the solu- 
tion of 11 linear equations of the type of equation (47)^ or (49) 

with one unknown and the interpolation of four points, which can be done 
in a very short time (of the order of 1 hour) when the net is drawn and 
the coefficients of the equations are determined. The pressure distri- 
bution obtained is presented in figure 7 £^d is conpared with the pres- 
sure distribution obtained from the exact method. 

For the back part of the body where the flow differs sli^tly from 
parallel flow, a cylindrical body with uniform flow at different Mach 
numbers can be considered as the basic body. In each region of the 
flow the Mach number for the basic flow is constant; however, the com- 
ponent u can be approximated conveniently by changing the basic- flow 
Mach number. 

When the body has a tail, a conical solution as proposed in refer- 
ence 6 for flow inside a tube can be assumed as the basic body. In this 
case, each streamline of the conical solution can be considered as the 
shape of the basic body (fig. 8), and the flow field can be obtained 
from conical- flow calculations, which can be determined easily in the 
hodograph plane. By changing the strength of the final shock of the 
conical solution, different ratios between maximum cross-sectional area 
and tall area can be obtained. Any part of the streamline can be 
assumed as the basic body shape. 

The present method does not require the existence of a linearized 
solution and, therefore, can be applied also at hi^ Mach muabers. 

This method permits obtaining the shape of the shock wave and taking 
into account entropy variations. Hi^ precision can be obtained by 
using several basic flow fields for the different regions of the body 
considered. Because of the slnplicity of the calculations, the system- 
atic calculations and tabulation of coefficients of equations ( 58 ) for 
different cones and different Mach numbers would be of great practical 
interest. 

Tabulated values can, be obtained also in the following way: For 

each cone OC of cone angle cp^ considered, a supeiposed flow must be 

calculated as shown, for exanple, in figure 9(a)- The values of 
and cp^ for this supeiposed flow are obtained at given points of the 

characteristic net. Because of the linearization of equations, if the 
superposed flow changes in intensity, the values of V-j^ and cp^ at 

every point change proportionately. Because of the conical property, 
if the point A (fig. 9(^)) nioves along OC, the flow field changes in 
scale; therefore, the effect on a point E due to the superposed flow 
field starting at A is equal to the effect of a linearized flow 
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field starting at A* and of Intensity at a corresponding 

point E* defined "by 


V ^ V ^ V 

^E ^A y-E 


Therefore, ■when the flow field for the disturbance AB is determined, 
the effect of any disturbance of the type of the disturbance AB in the 
entire flow field can be obtained. Then any body shape can be con- 
sidered as a supeiposition of flow fields of the type of flow corres- 
ponding to the disturbance placed along the basic cone. From -the 

sin5>le calculation of the flow for 'the shape AB, the velocity can be 
determined by means of "the equation 


V 


= Vo 


+ 




1 


'n 


■vdiere n is "the number of the superposed flows -that affect the point 
considered. 

Consider, for example, figure 10 . Several linearized flow dis- 
turbances must be superposed on -the basic conical flow field. First, 
a superposed conical flow at 0 that can be ob'tained from conical- flow 
calculations must be considered. At A^ a superposed disturbance must 

be added in order to satisfy the boundary conditions at A. If -the 
calculations have been performed for the disturbance at a^ a-nd 

for = 1 , the velocity at A can be ob'tained from 


Ta = (Vo)^ + (Vl)^ 


where (Vq)^ velocity of the basic flow at A, (Vp)^ is the 

velocity disturbance at A due to the conical flow superposed at 0 , 
and (^2)^ is proportional to 
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and can "be obtained from 


(V.) ^ 

^ (^*)a 



"Wiiere “b is determined "by 


y-b 


In a similar way determined at the points B, C, and so forth. 

At B another linearized flow field having velocity components 
and cp_ must he considered where 


('^3)2 ' - 


+ Wj * (^2)3 


and 


("3)b - 




where h is defined hy 

^h yh 


Because of the rapidity of calculation, the variation of any geometrical 
parameter can he investigated in practical applications without the 
necessity of a large amount of numerical work. 
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Conical Flow Field without Axial Symmetry 

The calculation of slender bodies without axial symmetry req.uires 
the determination of conical flow without axial symmetry, which can he 
done by means of the linearized characteristic method. The basic prob- 
lem of the determination of conical flow consists in determining the 
shape of the conical shock wave produced by the body. When the shape 
of the shock is determined, the flow field around the body can be 
obtained by meems of numerical calculations (see, for exan^ole, refer- 
ence 4 ) . Because the relation between the shape of the body and the 
shape of the shock wave is not known a priori, the method of linearized 
characteristics can be particularly useful for flow determination of 
this kind. An approximate shape of the shock wave is assumed as the 
basic solution and the flow field inside the shock is determined; then 
a linearized flow field is superposed in order to satisfy the boundary 
conditions at the body. The calculations are sinple if the basic flow 
can be determined analytically or numerically without a large amount of 
calculation. For exanple, for slender bodies the basic flow can be the 
axially symmetric flow for idilch values are available in tabulated form. 
Consider a conical shock wave which can be defined in polar coordinates 
as 


00 

’►s = (’I'o), + L 


00 

cos ne + ^ (ifj 
0 ® 


sin m 0 


(60) 


where all the values of and are small so that terms of the 

order of can be neglected. Such a shock wave is approximately of 
circular cross section, as is found for slender conical bodies. If the 
flow is assumed to have a symmetry plane, the second summation of equa- 
tion (60) is equal to zero. 


The velocity conponents in the radial direction in the 


tangential direction 


K)i> 


and in the direction normal to the 


shock shock wave are (see fig. 11); 


cos ilf 

(Vjj)^ = Yi sin Y cos a 
(vt)^ = sin Y sin a 


> (61) 


\ 
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■wiiere all the velocity ccm^jonents are referred to the limiting velocity 
and a is the angle of the dihedral between the plane nomnal to the 
shock wave and the plane containing the reference axis. Across the 
shock wave the following relations are valid: 




r - 1 

7 + 1 





7-1 
7 + 1 





- V^^sin^ sin 



- y - ^ (1 - Vi^cos^ - 

\ ny2 7 + 1 Vj sin i}f cos a 

(^t)i = (^ 1)2 
= (^r)2 ■ 


Vn ^sin^ 


sin^a) 


( 62 ) 


If the velocity con^onents v^,, v^^, and 

coordinates are considered at each point 
(fig. 11), 


w of the flow field in polar 


Tjf = + ^ n'ljr^ cos ne 




( 63 a) 


y -I (1 - Y-.^cob^'^ - V,^sln^^ sin^'j . 

K) , = - ^ ^ v. sin t sin^ 

\ “^Tjr 7 + 1 Vj_ sin f 1 


( 63 b) 
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(w)Tjr 


- sin ■'{r sin a cos a + 


7-1 
7 + 1 






sin 




sin a 

sin ilr cos a 


(63c) 


"but 


a 


dilf 

sin d0 



nijr^^ sin n0 
sin ■'}f 


Therefore^ the velocity con^jonents "behind the conical shock -wave of 
equation (60), if terms of the order of are neglected, are 


(^r)^ = Vi cos()Jro)^ - sln(tQ)^ ^ (t^)^ cos n0 


y - 1 1 “ ^l^cos^(^lr^)^ 

7+1 Vi sln(ir^) 


K) 




>(64) 


^cos(tJ |27i- 


1 - ^ 


Vi Bto (y„) 




cos n0 


7 - lA ■ 

7 + 1\ Vi sin^o„ 


(v)^s = 


Vi - 






sin n0 
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Then, th.e velocity con^jonents “behind the shock defined hy equation ( 60 ) 
can he expressed as: 


“rv 


' N),. H ne 


^ (65) 




•where “'^rQ^ Independent of 0, 

v_ and Vji correspond to the flow field for circular conical 
■'■o o 

shock , and ^v^.) , ("''’n) ■> circular cone ‘' 1 '' = ‘'I''o 


'n ' 'n 

are independent of n and are defined hy 


■ (''“oV 

I- -*'»'Os '*'Ob 


n>~i -vi , / i“^T ^cos^( 


°si _ 


i^^Os / ^ 




= Vi + 


n. 


^ \sin 




and 


. St 


'to. 


-V^ 


^ro + X '*'0 


1 - 




Jtn 
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It can te shawn that the flow field defined as in eguation (65) satis- 
fies the boundary at a surface of a conical body defined by 


= '•'oc ^ 


(67) 


if the terms are small and quantities of the order of or 

higher are negligible with respect to terms of the order of "' 1 ^^. At 

the surface of the body the velocity must be tangent to the surface of 
the body; therefore, at each point A of the body 


Y. 

sin de sin 

and, therefore, the boundary conditions can be expressed in the approxi- 
mation considered as 

— ^ 4 

sin tItq 



where the subscript ilr^ indicates q.uantities at the surface of the 

basic circular body end indicates quantities at the surface of the 
conical body considered where the parameters (^n)g given by equa- 
tion (60) and define the shape of the shock. Because each teim of the 
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ri^t-side summation is of the order of bou n dary conditions 

ar^ 


(’■no) 



= 0 







> 


(68) 


Equations ( 68 ) show that the basic flow defined by the components 

^o 

and is the flow corresponding to a circular cone of angle ijr^ 

and permits relation of the equation of the con i cal shock to the equa- 
tion of the conical body. The coefficients , ("%) ^ ('w)n 

for different values of n can he determined for a given free- stream 
Mach number and value of (^o)q (^o)g^ therefore, from equa- 

tions (60), (65), (67), and ( 68 ), the flow field around any conical body 
of the type given by equation (67) can be obtained idien the terms (^n)p 

are small. 


The determination of the quantities (Vp^ , ^Vj^^ , and ('w)^ as 

functions of ''I'’ can be obtained from the following equation (see refer- 
ence 3); 


M2- ^ 


v„2 + 




+ Vj, cot Tjr + ^1 - 


8 ilr \ a 2 / Bin ijr d© 


^^fi - 4 1 - 


8 v, 

j-L i yn ^ Ll 

g _2 sin ■'I'’ 80 


^ j. Tin 1 ^ 0 


(69) 
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Eq,uation. (69), because of expressions (65), gives 


N)n 




St 


V ‘ 

II - ^ 


-(^n)n 


/2v y „ 1 \^T + ■'^n cot t 

^ - ° o 

\^o ) y 2 

1 - 


y / V 1 + ■Vr, cot t \ 

(v 1 + 21^ J^o ^o : V V 2\ _ 

^ ^JrLl k 2 


n^„ 




o ^o sin t 


(70) 


Because of the approximation considered, the entropy remains constant 
in each meridian plane outside of a vorticose layer of infinitesimal 
thickness around the body (see reference 3). Indeed, 

^ H = sin t - v^ - v^ + v^v sin t + v^v cos t (71) 


and 


Vj, sin ijf — = - w — 

"" St S0 


(72) 


From equation (71) there results || of the order of (t^) j there- 

Is of the order of tj^^ and can be neglected. 


fore, -where Vj^ ^ 0, 
Then 


— — = 0 

rR St 


= -V, 


St 


-'w — + - + v^v„ - W^cot t ( 73 ) 


St St sintSe 
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or, in the approximation considered. 


_ 


m 


Sv-n 

Then is the radius of the hodograph diagram, in the plane 

0 = Constant (see reference 4) and and can he obtained 

from a step-hy-step calculation frcm ''1'' = '^''o to “ ^On means 
of the equations 




y (75) 


•where at station ilf is obtained from equation (70) and 


[(B)n]y = 






St 


(76) 


t 


can be calculated from the values of (^)n 


The value of Kw)Jl can be obtained from 

11 




( 77 ) 


and I can be obtained from equation (71) where 

\^t 
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a 


= So + Si ^ cos 


n 0 


( 78 ) 


and Si is independent of ijf and can "be deteimined from the equations 
of the shock from the expression 


S - §3- 




°s 


■where ''l^g is the inclination of the shock. Then 
-2 ^ 


7 R 


S. = T_ sin i fv-) + Vj, (vj,) + 

1 °-o ^\jf ^o' -^'n -^o' -^'n 


. ^r^^n sin Ilf + cos t 


( 79 ) 


while^ if the quantity Si is neglected and the flow is considered 
potential flow. 


Wn 


(^r)n 
sin ■')f 


( 80 ) 


The method presented has heen applied to the determination, for the con- 
dition of zero angle of attack, of the flow field around a cone ha-ving 
an elliptical cross section with axes in the ratio of about 1 to 3^ for 
•vdilch e3q)erimen‘tal da-fca ■were available at M = 1 . 8 , and for a cone of 
elliptical cross section of ratio 1 to 1 . 88 . The calculations have 
been peiformed in the following ■way: The bodies are shown in fig- 

ures 12(a) and 12(b). The value of i|r at 6 = 0 ° is equal to 6.3°, 
■while the ■value of ^If at 0 = 90° is equal to l 8 . 1 l-°. 

The angle ijr can be expressed as 


1)^ = cos 20 + ^ 2 ^ cos 40 + ■'If^ cos 60 + 


■'jf]^ cos 80 + if^^ cos 106 


(81) 
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i^2 


but 

y tan tn 
tan = tan f + 

COB% 

therefore, by determining the veO-ue of ijf at six points, the following 

values have been obtained: If tITq = 10°, 

C./ 

% " -5-07° 

Ilf = -0.84° 

3c 

\ = 0.43° 

. % = -0.16° 

5c 

Therefore, a, 10° cone at zero angle of attack is assumed as the basic 
body. Reference 5 gives tabulated values for a 10° cone at M-j_ = I.816 

and the calculations have been performed at this Mach number. The 
table gives 

iL = 34.45° 

°S 


The entropy variations S]_ are smal 1 and are neglected. Then, 

(^r)^ = 0.0255 
= 1-285 
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The values of , (v ) , and (w)„ hetweeii ilr and f have 

\ -^/n \ “ Og Oq 

been determined by means of equations (TO)^ ( 75 )^ and ( 80 ). Th^n ijr^ ^ 

S 

'^2q> ^ 5 g li 8 ’'ve been determined from equations ( 68 ) -where 




-2 X 0.6000 


The -ralues of ■'jfg obtained are: 


/ 


% = -0.24° 



5.1° X 10“3 


to = -1.7° X 10 -^ 
3s 

t, = 7.4° X 10 “^ 


% = -3° X 10" 


As is shown from the analysis of the -values of tq and tp, the shock 

is very close to a circular shock -wave even for large departure of the 
body from the circular cross section, and the effect of the terms corre- 
sponding to n = 6 , 8 , and 10 is very small. The velocity con 5 )onents 
at the surface of the body are obtained from equation ( 65 ) at t = 

and the pressure distribution presented in figure I 3 is ob-fcalned. In 
the same figure, the pressure distributions obtained by using -fche same 
calculated values of ("''n) ^ ^ around 

an ellipticad cone -with a cross section ha-vdng axes in -the ratio of 1 
to 1.88 are also shown. The conical body having an ellipse of axis 
ratio 1 to 3 has -the seime cross-sectional area as a circular cone 
of = 11°. Its pressure drag ob-tained from this calculation 

is Cjj = 0.099^ in con 5 )arison with 0.12 for the circular cone. The 

conical body having as cross section an ellipse -with 'axes in the ratio 
of 1 to 1.88 has a drag coefficient of O.IO 3 , -while the equi-valent 
circular cone of to = 10° 30* has Cj) = O.II 5 . Therefore, -those 
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calcula'tlons ludlca'be tlia'b conical 101168 of circular cross section 
have larger drag than cones of elliptical cross section. 

The results obtained, agree ■well ■with the experimental results, 
also, if the body shape chosen requires large values for the 
angle With the same flow fieMs ('^r)j^-> 

other conical shape having two planes of symmetry represented, by an 
equation (8I) "when ■'Jr^ is 10 ° can be obtained, at M = I.816. If the 

flow has only one plane of symmetry, only the terms in cos 0, cos 3 ^, 
Eind. so forth, must be consld.ered; wiiile, if no symmetry exists, terms 
in sin n0 and. cos n0 must be consid.ered. 

The flow fields defined by 

obtained a^nd given in tabulated form -without a large amount of numerical 
work as for circular cones, and, -therefore, the determination of conical 
bodies can be performed -wi-thout difficulty in a very short time. 

When the shock shape is some-what different from a cone having 
circular cross section, the basic flow field must be different from the 
axially symmetric. However, if the basic conical- flow components are 
expressed in the fonn 


v_ = v._ + v_ f(0) 


^ “b 


w = 


the basic flow can still be obtained by sol-ylng numerically the equa- 
tions of motion in two meridian planes, and, therefore, the basic flow 
can be determined exactly. For conical flow the linearized method can, 
then, have "wide application to any form of boundary conditions. 


Flow Fields Around Slender Bodies without Symmetry 

When the conical flow is determined, the method of characteristics 
can be applied to the determination of slender bodies. The equations 
used are similar to the equations for circular bodies at angles of 
attack, and can be directly derived from those equations (reference 2). 

Few values of n are required for the determination of the flow 
field, and one set of calculations can be used for several bodies having 
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the same haslc hody; therefore, the method can he of Interest for prac- 
tical applications. 


Quasi-Two-Dimensional Flow Fields 

In many general three-dimensional flow fields of practical Interest 
the flow is not too different from a two-dimensional flow, and, there- 
fore, the velocity field and entropy field can he expressed as in equa- 
tions (1) and (2) with good practical approximation. Flow fields of 
this kind are found, for exanq)le, in wings having plan forms which can 
he considered close to the two-dimensional type with some twist or a 
variation of thickness distrihutlon along the span. Flow fields of 
this kind can he considered also in some problems in which interference 
between a wing and a two-dimensional tall (downwash effects) or between 
a two-dimensional wing and a hody is considered. In all these problems 
of practical interest for the aiiplane design., the coiiqjonent w in the 
direction of the span of the wing can he considered small; therefore, 
equation (l) can he used and the con5>onents u^ and v^ depending on 

the three-dimensional effect can also he considered small. ' 

Equation (l?) expressed along the characteristic line of the first 
family = tan in the plane z = Constant is; 


Vo 


dYi 

^ . tan Po 


dx 


sin^PQ dS]^ 
7R . dx 


9^Bi 


+ IlC -'^^11 '*^ PQ PQ 

Vq 1 dz Vo cos(cpp + P^) 


(82) 


while along the second characteristic line X2 = “ ^o) 

following equation is valid; 


Vq dx 


+ tan Pq 


^1 

dx 


sin^o dSi 
7R dx 


- 9^32 + 


Vo 


C2 = 


d^ tan Pq sin p^ 
iz Vq cos(cpQ - Pq) 


(83) 


where C^, and are defined by equations (20). Along each 

streamline Sq, 

dw^ StVi a^2 

dsQ dz dz 7 RVq 


m 
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Equations (82), ( 83 ) ^ and (8^)-) permit the determination of the flow 
field "by relatively single procedui^s. 

Consider, for example, a wing having twist, variable profile dis- 
tribution, and variable chord, as shown in figure l4. The wing can be 
analyzed by means of the linearized characteristics method in the fol- 
lowing way; First, the root and tip profiles are considered. Section a 
and section b have different relative thicknesses and chords. 

If the variation from a to b is linear, the properties of a 
two-dimensional cross section at any station c can be obtained by 
means of linear interpolation between the corresponding values at a 
emd b. 


The profiles a and b are analyzed by means of two-dimensional- 
flow theory and the characteristic net, and the values of the coef- 
ficients B 2 , and C 2 are determined from two-dimensional 

considerations. If entropy effects are neglected or Incorporated in 
the linearized flow, the coefficients B 2 , and Cg can be 

determined as for the case of two-dimensional potential flow at each 
point of the axis (equations ( 29 ), ( 30 ), and (3^)) and are constant 
along characteristic lines of the first family. Then the linearized 
flow is defined as the flow that considers the three-dimensional effects 
and the entropy distribution. Therefore, 


u = Uo(x,y,z) + aui(x,y,z) 
V = Vg(x,y,z) + av^(x,y,z) 
w = aw^(x,y,z) 


( 85 ) 


where Uq and Vq are the potential- flow solutions in the pleine 

z = Constant arfl satisfy the boundary conditions in the plane 
z = Constant, Uj_ and Vq^ are the conq)onents due to the presence 

of w^^ (and of the variation of entropy), emd a cem be a coefficient, 

for exan 5 )le, proportional to the twist distribution or to the thickness 
variation. Because Uq and Vq are functions of x, y, and z, while 

for the basic flow they have been determined from two-dimensional con- 
siderations, and are not zero; therefore, equation (22) 

oz oz 

8w 3wu,8wv u .TtSY 

8sq 8xT 8yY Y * Y 8z 


becomes 
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or 


= ^ + 
SSq Sz Sz 


+ 


Sz 



7 RV 0 


( 86 ) 


ov 

vhere — - is the variation of the velocity conrponent for the basic 
dz 

flow. 

T'Jhen is considered small, in all the flow field the terms 

2uw ^ and — can still he neglected in the differential equa- 

a2 Bz a2 Sz 

tions along the characteristic lines, and, therefore, equations (82) 
and ( 83 ) are still valid. 

At each plane z = Constant, the characteristic net is kaown; 
therefore, the intersections of the shock wave for the total flow with 
a plane y = Constant can he determined from characteristic calcu- 
lations. If 0 is a point at the leading edge of the wing (fig. llf(h)), 
the shock wave at 0 can he obtained from shock-wave considerations and 
from the boundary conditions because 5 at 0 is known, and at 0 the 
shock is two-dimensional. Therefore, the velocity conponents U 3 _ 

and V-, at 0 are zero, while is given by equation (86). If the 

OSq 

plane x = is assumed to be close to the plane x = x^, the charac- 

teristic lines BA. and CA can be drawn in any meridian plane considered 
for the basic flow. 

At the point 0, Vp is zero €tnd, in the nel^borhood of 0 along 
the shock L the velocity can be expressed as 

(Yl)^ - AL (87) 


Now, along the shock wave the direction of the velocity behind the shock 
is related to the intensity from the equations of the shock wave; there- 

SCDn 

fore, the direction cp.^ = along the shock and the shape of the 

1 oL 


8Vn 8Vn 

shock as a function of -— =■ are determined when 

oL oL 


along the shock 


is known. 
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If the velocity at a point A of the hody is 




As 


(88) 


then the velocity at C and B can he determined as a function 

of — from equations ( 82 ) and (83) ^ because the value of w-, at A 
os 


is given from equation ( 86 ) and is knovm., and the value of 
obtained from the value of -w-j^ at A in several planes z 


, 

— can be 
dz 

= Constant. 


If in equations ( 82 ) and (83) the values of and cp^ at 

B and C are expressed by means of equation (87) and of the equations 
of the shock waves -vdiich give the coefficients of the expressions 



SL dx 


8<Pi 8cp2^ 

8x 8x 


then 

and. 


equations 

therefore. 


(82) and (83) give two relations between 
8Vi SVi 

8s 8 l 


can be determined. 


8s 


and 



The equation of the shock wave can relate u as a function of v, 
or as a function of cp^ i-rithout the necessity of the component w, 

because the con^jonent w is proportional to the inclination t] of the 
tangent to the shock with the plane x = Constant, eind and cp^ are 

functions of cos tj, but cos tj = M^l - ~ ^1 ^ approxi- 

mation considered here (fig. 14 (c)). The con5)onents u and v at B 
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the position of B having heen determined, in each merid i an plane 
z = Constant, the intersection of the shock in the plane y = Constant 
is obtained and w at B determihed. Then a point D is interpolated 
in each meridian plane and the point E is obtained. Then point E is 
determined. In order to obtained ('w)ji, the stre am ! 3 , n e DF’ for the 

basic flow must be drawn and interpolated between F* and E. 

\Sz /o 

In a similar way, all the flow field can be obtained. The line TT* 
defining the plan form in figure l4(a) must be outside of the Mach 
conoid from T. Because few points along each profile are required, the 
largest amount of work for such a calculation is represented by the con- 
struction of a basic characteristic net which permits obtaining points 
which sin^jlify the determination of w. By ch a nging the value of the 
coefficient a, different thickness distributions or different twists 
can be considered. The new distributions must be obtained by changing 
proportionately the variation of thickness or twist with respect to the 
basic wing and by varying in proportion the value of a. 


CONCLUDING EMABKS 


The method of characteristics for supersonic flow has been sin^)!!- 
fied by assuming that one of the velocity con5)onent8 or the effect on 
the velocity conponents due to variation of one physical parameter is 
small, so that the square of the velocity ccanponents considered small 
can be neglected. By means of this simplification, the flow field can 
be represented as the supeiposition on a basic flow field (which is not 
linear and must be deteimJ.ned by the method of characteristics) of 
linearized flow fields which are defined by a differential equation 
with variable, but known, coefficients. 

The calculations of these linearized flow fields can be perfonmed 
along the chaiacterlstlc net of the basic flow field. The method has 
been applied (a) to the two-dimensional flow with entropy gradient, 
which has been transformed to a basic potential flow on which a line- 
earlzed flow due to the entropy gradient is superposed, (b) to axially 
symmetric problems where conical or cylindrical flows are considered 
as the basic flow, (c) to the determination of the flow field around 
cones or slender bodies without axial symmetry, and (d) to particular 
three-dimensional flows which can be simulated as a basic two- 
dimensional flow on Tdiich three-dimensional linearized flows are super- 
posed. Application (b) permits obtaining in a simple way the flow 
field aroimd bodies of revolution without using linearized theory, and 
indicates the possibility of using tabulated values for such determi- 
nation. Application (c) penults the determination of flow fields not 
yet determined by the method of characteristics. Any such conical 
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flovra can "be determined ’by using tabulated values tlxat can be obtained 
as for cones of circular cross section at small angles of attack. The 
application in (d) can be of interest for ■wings of approximately -two- 
dimensional form ha-ving -twist or -thickness -mrlatlon along -fche span and 
to interference problems. 


Langley Aeronautical Laboratory 

National Ad-visory Committee for Aeronautics 
Langley Field, Va., July 24, 1951 
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